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We report a Dalitz plot analysis of charmless hadronic decays of charged B mesons to the final state 
AgTr^TT® using the full BABAR dataset of 470.9 ± 2.8 million BB events collected at the T(45') reso¬ 
nance. We measure the overall branching fraction and CP asymmetry to be B ^ = 

(45.9 ± 2.6 ± 3.0t^;g) x 10“® and Acp {B+ K°TT+n°) = 0.07 ± 0.05 ± O.OSlK g^, where the un¬ 
certainties are statistical, systematic, and due to the signal model, respectively. This is the first 
measurement of the branching fraction for S'* ** " —>■ We find first evidence of a CP asym¬ 
metry in B+ A*(892) + 7r® decays: Acp {B+ A*(892)+7r°) = -0.52 ± 0.14 ± 0.04t° ®^. The 

significance of this asymmetry, including systematic and model uncertainties, is 3.4 standard devi¬ 
ations. We also measure the branching fractions and CP asymmetries for three other intermediate 
decay modes. 

PACS numbers: 13.25.Hw, 12.15.Hh, ll.30.Er 


I. INTRODUCTION 

The Gabibbo-Kobayashi-Maskawa (CKM) mecha¬ 
nism [Ei for quark mixing describes all weak charged 
current transitions between quarks in terms of a unitarity 
matrix with four parameters: three rotation angles and 
an irreducible phase. The unitarity of the CKM matrix 
is usually expressed as triangle relationships among its 
elements. The interference between tree-level and loop 
(“penguin”) amplitudes can give rise to direct CP vio¬ 
lation, which is sensitive to the angles of the Unitarity 
Triangle, denoted a, /3, and 7 . Measurements of the pa¬ 
rameters of the CKM matrix provide an important test 
of the Standard Model (SM) since any deviation from 
unitarity or discrepancies between measurements of the 
same parameter in different decay processes would im¬ 
ply a possible signature of new physics. Tree ampli¬ 
tudes in B ^ K*tt decays are sensitive to 7 , which can 
be extracted from interferences between the intermedi¬ 
ate states that populate the ATtttt Dalitz plane. How¬ 
ever, these amplitudes are Cabibbo-suppressed relative 
to contributions carrying a different phase and involv¬ 
ing radiation of either a gluon (QCD penguin) or photon 
(electroweak penguin or EWP) from a loop. 
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Perugia, Italy 

^Now at: Laboratoire de Physique Nucleaire et de Hautes Energies, 
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QCD penguin contributions can be eliminated by con¬ 
structing a linear combination of the weak decay am- 
plitudes for 5+ —^ K*Tr to form a pure isospin / = § 
state 0 : 

Ab = A{K*°Tr+)+V2A{K*+Tr°) . (1) 

Since all transitions from / = ^ to I = | states occur via 
only A/ = 1 operators. As is free from QCD contribu¬ 
tions. The weak phase of As is often denoted as 

2 

$3 =-lArg(A3/A|), ( 2 ) 

where As is the CP conjugate of the amplitude in Eq. O- 
The phase <i)| in Eq. © is the CKM angle 7 in the 
absence of EWP contributions 0 . 

Measurements of the rates and CP asymmetries in 
B ^ Ktt have generated considerable interest because 
of possible hints of new-physics contributions @,0. Of 
particular interest is the difference, AAqp, between the 
CP asymmetry in 5+ —>■ iA+ 7 r° and the CP asymme¬ 
try in ^ Ar+ 7 r“, which in the SM is expected to be 
consistent with zero within the theoretical uncertainties 
assuming U-spin symmetry and in the absence of color- 
suppressed tree and electroweak amplitudes 0,0. Using 
the average values of Acp of AT+tt^ and A+tt” decays 0 , 
Is.A cp (Ktt) is 

AAcp {Ktt) = Acp — Acp (A'''7r“) 

= 0.122 ± 0.022, (3) 

which differs from zero by 5.5 standard deviations. Un¬ 
fortunately, hadronic uncertainties prevent a clear inter¬ 
pretation of these results in terms of the new-physics 
implications 0113 . Additional information can be ob¬ 
tained through studies of the related vector-pseudoscalar 
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decays B K*Tr and B —7> Kp [iMl , for which the ra¬ 
tios of tree-to-penguin amplitudes are expected to be two 
to three times larger than ior B ^ Ktt decays. Hence, 
B —)> K*Tr and B ^ Kp decays could have considerably 
larger CP asymmetries. 

In this article, we present the results from an ampli¬ 
tude analysis of —>■ decays. The inclusion of 

charge conjugate processes is implied throughout this ar¬ 
ticle, except when referring to CP asymmetries. This is 
the first Dalitz plot analysis of this decay by BaBAR; the 
only previous BaBAR analysis of this decay was restricted 
to measuring the branching fraction and CP asymmetry 
of B^ —7> K^p'^ [I^. An upper limit on the branching 
fraction for -A was set by the CLEO Col¬ 
laboration: B{B^ -A < 66 X 10“® [l^ . 

Two contributions to the Ar°7r+7r° final state arise 
from the resonant decays B^ -A A'*(892)°7r+ and 
B+ -A A:*(892)+7r°. Although both the rate and CP 
asymmetries for A'*(892)°7r+ have been well mea¬ 

sured, with Ar*° —AT+tt”, by both the BaBAR [l^ and 
Belle [13 Collaborations, the measurements of the rate 
and CP asymmetry for 5+ —Ar*(892)+7r° [l^ have sig- 
nihcant statistical uncertainties and could benefit from 
the additional information provided by a full amplitude 
analysis. In Table U we review the existing measurements 
of the rates and CP asymmetries in the B -A Ar*(892)7r 
system. 


TABLE I: Average values of the branching fractions B and 
CP asymmetries Acp for B -A A*(892)7r decays as deter¬ 
mined by the Heavy Flavor Averaging Group [^. 


Mode 

5(10"'’) 

Acp 

References 

A*+7r“ 

8.5 ±0.7 

-0.23 ±0.06 

[19-22] 

A*+7r“ 

8.2 ± 1.8 

-0.06 ±0.24 

[18] 

K*°tv+ 

Q Q +0-8 
y-y _o.9 

-0.038 ± 0.042 

[16, Uj 

O 

O 

2.5 ±0.6 

-0.15 ±0.13 

[19, 2^ 


This article is organised as follows. The isobar model 
used to parameterize the complex amplitudes describing 
the intermediate resonances contributing to the AT^Tr+Tr^ 
final state is presented in Section El A brief descrip¬ 
tion of the BaBAR detector and the dataset is given in 
Section EH The event reconstruction and selection are 
discussed in detail in Section llvl the background study 
in Section El and a description of the extended maxi¬ 
mum likelihood fit in Section lyll The results are given in 
Section [Vlli and a study of the systematic uncertainties 
is presented in Section IVIIII In Section IIXI we provide 
a summary and conclusion, discussing the results and 
combining the branching fractions and CP asymmetries 
for the decays B+ -A K*{892)°tt+, B+ -a K^{U30)°tt+, 
and Ar*(892)+7r° with previous BaBAR results 

obtained from the final states Ar+7r“7r+ and 

B+ -A AT+ttOttO. 


II. AMPLITUDE ANALYSIS FORMALISM 


A number of intermediate states contribute to the 
decay -a AT^Tr+Tr^. Their individual contributions 
are measured by performing a maximum likelihood fit 
to the distribution of events in the Dalitz plot formed 
from the two variables, mCo + and n- We use the 

Laura++ [2^ software to perform this fit. 

The total signal amplitudes for the B~^ and the B~ 
decays are given in the isobar formalism by [2^, [2^ 

A (^KO^+’^i+^o) {^KOn+’^l+no) . (4) 

3 

A ’ ™'7r-7i-o) (jTLkOtt- ’ (®) 

3 

where Cj is the complex coefficient for a given resonant 
decay mode j contributing to the Dalitz plot. This com¬ 
plex coefficient contains the weak-interaction phase de¬ 
pendence that is measured relative to one of the con¬ 
tributing resonant channels. In this article we report 
results for the relative phases between each pair of am¬ 
plitudes. 

The function Fj describes the dynamics of the decay 
amplitudes and is the product of a resonant lineshape 
(i?j), two Blatt-Weisskopf barrier factors (Xl), and 
an angular-dependent term 

Fj = Rj X Xl{\p\ , IpoI) X XlUcH , |go|) x , (6) 

where L is the orbital angular momentum between the in¬ 
termediate resonance and the bachelor particle (the bach¬ 
elor particle is the daughter of the B decay that does not 
arise from the resonance), q is the momentum of one of 
the daughters of the resonance in the rest frame of the 
resonance, p is the momentum of the bachelor particle in 
the rest frame of the resonance, and po and qo are the 
values of p and q, respectively, at the nominal mass of 
the resonance. The Blatt-Weisskopf barrier factors are 
given by 


u , uo ) = 

1, 

(7) 

Al=i(|u|,|uo|) = 

/l ± Zq 

y 1 + z’ 

(8) 

Xl=2{\u\ ,\uo\) = 

1 (zo - 3)^ ± 9zo 
y (z — 3)^ ± 9z 

(9) 


where z = (|m| rew)^; = (If^ol ^bw)^, u is either q orp, 
and tbw = 4.0 (GeV/c)~ is the meson radius parameter. 
The uncertainty in rewj used for systematic variations, 
is ±2 ( GeV/c)~^ for the K* resonances, and ranges from 
-1.0 to -f2.0 (GeV/c)”^ for the p(770)+ [2|. The an¬ 
gular term depends on the spin of the resonance and is 
given by m 113 




= 1 , 


-2p.q, 

4 r 9 
o 3(p.<^ 


(Ip1 \q\f 


( 10 ) 

( 11 ) 


( 12 ) 
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(Kg7r*);r° (Kg7r°)7i^ (7i^7r°)Kg 



FIG. 1: Schematic representation of the definitions of q and p 
used in this analysis for the (left) KsTt'^, (center) , and 

(right) 7r+7r° resonances. 


TABLE II: Parameters of the Dalitz plot model for 
B+ K%-k+'k° used in the nominal fit. The mass and width 
of the p(770)''' and their uncertainties are taken from the anal¬ 
yses by the ALEPH [3l| and CMD2 Collaborations. All 
other parameters are taken from Ref. The resonance 

shapes are a Gounaris-Sakurai (GS) function, a relativistic 
Breit-Wigner (RBW) function, or based on measurements by 
the LASS Gollaboration [s^, with a the scattering length and 
r the effective range of the LASS parametrization. 


Resonance Lineshape 

Parameters 

Resonance mass Width 
(MeV/c^) (MeV) 

P(770)+ 

GS 

775.5 ± 0.6 

148.2 ±0.8 

A* (892)+ 

RBW 

891.7 ±0.3 

50.8 ±0.9 

A* (892)° 

RBW 

896.1 ±0.2 

50.7 ±0.6 

(A^);°/+ 

LASS 

1412 ± 50 

294 ± 80 



mcutofi = 1800 MeV/c^ [1^ 



a = 2.1 ±0.1 (GeV/c)-° [16] 



r = 3.3 ±0.3 (GeV/c)-i [1^ 


The choice of which resonance daughter is defined to 
carry the momentum g is a matter of convention. How¬ 
ever, its definition is important when comparing mea¬ 
surements from different experiments. In Fig. [U we il¬ 
lustrate the momentum definitions used for the 
and tt+tt^ resonance combinations. 

Table H lists the resonances used to model the signal. 
We determine a nominal model from data by studying 
changes in the log likelihood values for the best fit when 
omitting or adding a resonance to the fit model, as de¬ 
scribed in Section ED 

For the AT*(892)° and A'*(892)+ resonances, we use a 
relativistic Breit-Wigner (RBW) lineshape 28|: 


E 


RBW 


(m) = 


1 


Trig — m? — imQV{m) ’ 


(13) 


where m is the two-body invariant mass and F(to) is 
the mass-dependent width. In general, for a resonance 


decaying to spin-0 particles, F(m) can be expressed as 

(lii) (14) 

where mg and Fo are the nominal mass and width of the 
resonance. 

The Gounaris-Sakurai (GS) parametrization [s^ is 
used to describe the lineshape of the p resonance decaying 
into two pions. The parametrization takes the form 

nGS ^ 1-bFo-d/mo 

^ rrig — vn? + f{m) — iTOor(m) ’ 

where F(m) is given by Eq. (fHl) . Expressions for f{m), 
in terms of Fq and m, and the constant d can be found 
in Ref. [^ . The parameters specifying the p lineshape 
are taken from Refs. [ID, which provides lineshape 
information derived from fits to e“'"e“ annihilation and r 
lepton decay data. 

For the component of the Ktt spectrum, de- 

noted we make use of the LASS parametriza¬ 

tion [ 3 ^, which consists of a Kg resonant term together 
with an effective-range, nonresonant component to de¬ 
scribe the slowly increasing phase as a function of the 
ATtt mass: 


nLASS _ 
Kj 


1^ cot 5 b - i \q\ 


d22(5s 


moToj^ 
^ ^\qo\ 


(mg-m2)-Imoro^’ 


(16) 


where cot5_B = ^ The values used for the 

scattering length a and the effective range r are given 
in Table HD The effective-range component has a cutoff 
imposed at ISOOMeV/c^ [l^. Integrating separately the 
resonant term, the effective-range term, and the coher¬ 
ent sum, we find that the Arg(1430)° and the ArQ(1430)''' 
resonances account for 88% of the sum, and the effec¬ 
tive range component 49%; the 37% excess is due to de¬ 
structive interference between the two terms. The LASS 
parametrization is the least-well-determined component 
of the signal model; we discuss the impact of these un¬ 
certainties in Section EUD 

The complex coefficients Cj and cj in Eqs. (I4l5l) 
can be parametrized in different ways; we follow the 
parametrization used in Ref. [TDj i as it avoids a bias in 
the measurement of amplitudes and phases when the res¬ 
onant components have small magnitudes: 


Cj = {xj+Axj) + i{yj+Ayj), (17) 

Cj = {xj - Axj) + i {yj - Ayj) , 

where Xj ± Axj and yj ± Ayj are the real and imagi¬ 
nary parts of the amplitudes. The quantities Axj and 
Ayj parametrize the CP violation in the decay. The CP 


















asymmetry for a given intermediate state is given by 



2 (xjAxj + VjAyj) 

x'j + Ax^ + yj + ' 


(18) 

(19) 


The results quoted for the resonances in the follow¬ 
ing analysis use fit fractions (FF^) as phase-convention- 
independent quantities representing the fractional rate of 
each contribution in the Dalitz plot. The FF for mode j 
is defined as 


FF . =__ I _ 

f f (l^l^ + l^r) dmj^^dml^ 


( 20 ) 


The sum of all the fit fractions does not necessarily yield 
unity due to constructive and destructive interference, as 
quantified by the interference fit fractions given by |30| 


pp.. ^ / f 2Re [c^c*F,F*] dm^^dml^ 

I I\J2k<^kFkf dmj^^dml^ 

The parameters xj, Axj, yj, and Ayj are determined 
in the fit, except for the reference amplitude. Fit frac¬ 
tions, relative phases, and asymmetries are derived from 
the fit parameters and their statistical uncertainties de¬ 
termined from pseudo experiments generated from the fit 
results. 


III. THE BABAR DETECTOR AND MC 
SIMULATION 

The data used in the analysis were collected with 
the BaBAR detector at the PEP-II asymmetric-energy 
e+e” collider at SLAG National Accelerator Laboratory. 
The sample consists of 429 fb~^ of integrated luminosity 
recorded at the T{iS) resonance mass (“on-peak”) and 
45 fb“^ collected 40 MeV below the resonance mass (“off- 
peak” ) [s^ . The on-peak sample corresponds to the full 
BaBAR T(4:S) dataset and contains 470.9 ± 2.8 million 
BB events [^. A detailed description of the BaBAR de¬ 
tector is given in Refs. [H, Charged-particle tracks 
are measured by means of a five-layer double-sided sil¬ 
icon vertex tracker (SVT) and a 40-layer drift chamber 
(DCH), both positioned within a solenoid that provides 
a 1.5 T magnetic field. Charged-particle identification 
is achieved by combining the information from a ring¬ 
imaging Cherenkov detector (DIRC) and specific ioniza¬ 
tion energy loss (dE/dx) measurements from the DCH 
and SVT. Photons are detected and their energies mea¬ 
sured in a CsI(Tl) electromagnetic calorimeter (EMC). 
Muon candidates are identified in the instrumented flux 
return of the solenoid. 

We use GEANT4-based software to simulate the de¬ 
tector response and account for the varying beam and 


experimental conditions [s^ [s^. The EvtGen and 
JetsetT .4 [4l| software packages are used to generate sig¬ 
nal and background Monte-Carlo (MC) event samples in 
order to determine efficiencies and evaluate background 
contributions for different selection criteria. 


IV. EVENT SELECTION 


We reconstruct B^ -A candidates from one 7r° 

candidate, one candidate reconstructed from a pair of 
oppositely charged pions, and a charged pion candidate. 
The TT^ candidate is formed from a pair of neutral en¬ 
ergy clusters in the EMC with laboratory energies above 
0.05 GeV and lateral moments between 0.01 and 0.6. 
We require the invariant mass of the reconstructed to 
lie in the range 0.11 < < 0.16 GeV/c^. The ATg can¬ 

didate is required to have a tt+tt” invariant mass within 
15 MeV/c^ of the mass [28|, and a proper decay time 


greater than 0.5 x 10“^^s. To reduce combinatorial back¬ 
ground, we also require that the ATg candidates have a 
vertex probability greater than 10“® and that the cosine 
of the angle between the ATg momentum direction and 
the ATg flight direction (as determined by the interaction 
point and the ATg vertex) be greater than 0.995. For the 
7r+ candidate, we use information from the tracking sys¬ 
tems, the EMC, and the DIRC to select a charged track 
consistent with the pion hypothesis. We constrain the 
7r+ track and ATg candidate to originate from a common 
vertex. 

Signal events that are misreconstructed with the de¬ 
cay products of one or more daughters completely or 
partially exchanged with other particles in the rest of 
the event have degraded kinematic resolution. We re¬ 
fer to these as “self-cross-feed” (SCF) events. This mis- 
reconstruction has a strong dependence on the energy 
of the particles concerned and is more frequent for low- 
energy particles, i.e., for decays in the corners of the 
Dalitz plot. Because of the presence of a 7r° in the final 
state, there is a significant probability for signal events to 
be misreconstructed due to low-energy photons from the 
7r° decay. Using a classification based on MC informa¬ 
tion, we find that in simulated events the SCF fraction 
depends strongly on the resonant substructure of the sig¬ 
nal and ranges from 34% for B+ -A A:*(892)+7r0 to 50% 
for B~^ -A p(770)+Arg. In events simulated uniformly in 
phase space, hereafter referred to as nonresonant MC, 
the SCF fraction varies from less than 10% in the center 
of the Dalitz plot to almost 70% in the two corners of 
the Dalitz plot, where either the 7r° or the 7r+ has low 
energy. We describe how the SCF events are handled in 
Section IVII 

In order to suppress the dominant background, due to 
continuum e“''e“ -A- qq {q = u, d, s, c) events, we employ 
a boosted decision tree (BDT) algorithm that combines 
four variables commonly used to discriminate jet-like qq 
events from the more spherical BB events in the e“'"e“ 
center-of-mass (CM) frame. The first of these is the ratio 
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of the second-to-zeroth order momentum-weighted Leg¬ 
endre polynomial moments, 

E ~ 1) Pi 

^2 zgROE _ (00\ 

To- ^ ^ 

iGROE 

where the summations are over all tracks and neutral 
clusters in the event, excluding those that form the B 
candidate (the “rest of the event” or ROE); pi is the par¬ 
ticle momentum, and 9i is the angle between the particle 
and the thrust axis of the B candidate, hereafter also 
referred to as the B. The three other variables entering 
the BDT are the absolute value of the cosine of the an¬ 
gle between the B direction and the collision axis, the 
zeroth-order momentum-weighted Legendre polynomial 
moment, and the absolute value of the output of another 
BDT used for “flavor tagging”, i.e., for distinguishing B 
from B decays using inclusive properties of the decay of 
the other B meson in the T(45') —^ BB event [d^. The 
momentum-weighted Legendre polynomial moments and 
the cosine of the angle between the B direction and the 
beam axis are calculated in the e+e“ CM frame. The 
BDT is trained on a sample of signal MC events and off- 
peak data. We apply a loose criterion on the BDT output 
of BDTout > 0.06, which retains approximately 70% of 
the signal while rejecting 92% of the qq background. 

In addition to BDTout , we use two kinematic variables 
to distinguish the signal from the background: 

’^ES = (23) 

AE = E*B-^/s/2, (24) 

where 

Ey^ = (s/2 -I- Pe+e- • Pb) jE^+e- , (25) 

and where ^/s is the total e“'"e“ CM energy, with 
(Ee+e-,Pe+e-) ^nd {Eb,Pb) the four-momenta of the 
initial e“'"e“ system and the B candidate, respectively, 
both measured in the lab frame, while the star indicates 
the e+e" CM frame. The signal toes distribution for cor¬ 
rectly reconstructed events is approximately independent 
of their position in the RlgTr+Tr^ Dalitz plot and peaks 
near the B mass with a resolution of about 3.4MeV/c^. 

We retain all candidates satisfying the following selec¬ 
tion criteria: 5.23 < toes < 5.29GeV/c^ and —0.3 < 
AE < 0.3 GeV. The signal region, where the final fit 
to data is performed, is defined by the tighter crite¬ 
ria 5.260 < toes < 5.287 GeV/c2 and -0.20 < AE < 
0.15 GeV. We also use candidates in the sideband re¬ 
gion of Toes defined by 5.23 < toes < 5.26GeV/c^ and 
—0.20 < AE < 0.15 GeV and subtract from distribu¬ 
tions for these events the BB background contributions 
predicted by MC simulations. We then add these dis¬ 
tributions to the off-peak data distributions to increase 
the statistical precision of our model of the Dalitz plot 
distribution for continuum background. 


Each of the B candidates is refit to determine the 
Dalitz plot variables. In these fits the iLg7r+7r° invari¬ 
ant mass is constrained to the world average value of the 
B mass [2^ to improve position resolution within the 
Dalitz plot. 

We find that 20% of the remaining events in nonreso¬ 
nant MC have two or more candidates. We choose the 
best candidate in multiple-candidate events based on the 
highest R-vertex probability. This procedure is found to 
select a correctly reconstructed candidate more than 60% 
of the time and does not bias the fit variables. 

The reconstruction efficiency over the Dalitz plot is 
modeled using a two-dimensional (2D) binned distri¬ 
bution based on a generated sample of approximately 
2 X 10® simulated —>■ MC events, where the 

events uniformly populate phase space. All selection cri¬ 
teria are applied except for those corresponding to a ATtt 
invariant-mass veto described below, which is taken into 
account separately. The 2D histogram of reconstructed 
MC events is then divided by the 2D histogram of the 
generated MC events. In order to expand regions of phase 
space with large efficiency variations, the Dalitz plot vari¬ 
ables are transformed into “square Dalitz plot” [44| coor¬ 
dinates. We obtain an average efficiency, for nonresonant 
MC events, of approximately 15%. In the likelihood fit 
we use an event-by-event efficiency that depends on the 
Dalitz plot position. 

V. BB BACKGROUNDS 

In addition to continuum events, background arises 
from non-signal BB events. A major source of BB 
background arises from R+ —>• D® (—>• AlgTr®) 7r+ de¬ 
cays. To suppress this background, we veto events with 
1.804 < mjio^o < 1.924 GeV/c2. 

The remaining BB backgrounds are studied using MC 
simulations and classified based on the shape of the toes, 
AE, and Dalitz plot distributions. We identify nine cat¬ 
egories oi BB backgrounds: categories 1, 2 and 3 include 
different types of three- and four-body B decays involv¬ 
ing an intermediate D meson; categories 4 and 5 include 
charmless four-body B decays to intermediate resonances 
where a tt® in the final state is not reconstructed; cate¬ 
gories 6 and 7 include two-body B decays with a radi¬ 
ated photon misreconstructed as a tt® decay product or 
where the tt® arises from the other B decay; category 8 
includes charmless three-body B decays where a charged 
pion is interchanged with a tt® meson from the other B\ 
and finally category 9 includes all other simulated BB 
background contributions. Within each category, each 
of the Toes, AE, BDTout, and Dalitz plot distributions 
are formed by combining the contributions of all decay 
modes in the category. The combinations are done by 
normalizing the distributions for each decay mode to the 
expected number of events in the recorded data sam¬ 
ple, which is estimated using reconstruction efficiencies 
determined from MC, the number of BB pairs in the 




10 


recorded data sample, and the branching fractions listed 
in Refs. For each category, the histograms of mssj 

AE, BDTout, and the Dalitz plot variables are used as 
the probability density functions (PDF) in the likelihood 
fit to data to model the BB background. 


VI. THE MAXIMUM LIKELIHOOD FIT 


The extended likelihood function is given by 


£ = exp 



(26) 


n 


'^NkVl AE,BDT out 1 


where Nk is the number of candidates in each signal or 
background category k, Ng is the total number of events 
in the data sample, and VI (the PDF for category k and 
event i) is the product of the PDFs describing the Dalitz 
plot, mESj and BDTout distributions, with qb the 
charge of the B candidate. 

To avoid possible biases in the determination of the 
fit parameters [45|, we use MC samples to study corre¬ 
lations between the fit variables and the Dalitz plot pa¬ 
rameters, TO^o^+ and We find that for correctly 

reconstructed signal candidates, the AE distribution is 
strongly dependent on 771 ^ 0 ,^.+ . This is mostly due to 
a dependence of the energy resolution of the B candi¬ 
date on the 7r° momentum. For SCF signal candidates, 
both the Toes and AE distributions depend on all three 
two-body invariant masses: rrixo^-K+i and TOtt+tto. 

The mESj AE^ and BDTout distributions for continuum 
and BB backgrounds have negligible correlations with 
the Dalitz plot parameters. 

For correctly reconstructed signal candidates, the toes 
and AE PDFs are parameterized by a Cruijff function, 
which is given by (omitting normalization factor) 


/cruijff(a;) = exp 


— (x — rnf' 


+ ^l,r (a; - to) 


(27) 


where to gives the peak of the distribution and the asym¬ 
metric width of the distribution is given hy a l for x < m 
and CTfl for x > m. The asymmetric modulation is sim¬ 
ilarly given by aL for a; < TO and or for x > m. The 
AE PDF parameters are calculated on an event-by-event 
basis in terms of the invariant mass, as a linear 

function for m‘^o^+ < 20GeV^/c'* and as a quadratic 

function for m‘^o^+ > 20GeV^/c'*. These functions are 
determined by fitting the AE distribution in large non¬ 
resonant MG samples. For the SCF signal, in order to 
follow the rapid shape variations across the Dalitz plot of 
the Toes and AE distributions, we divide the Dalitz plot 
into several regions as illustrated in Fig.[5J Each letter in¬ 
dicates whether the dependence is on to^+^o, or 

TO^o^o • The regions are chosen based on the distribution 


CM 

cT" 

o 25 


> 



FIG. 2: Diagram illustrating the division of the Dalitz plot 
into different regions for the definition of the PDFs for self¬ 
crossfeed signal events. Each letter indicates whether the de¬ 
pendence is on (A), rn']^o^^+ (B), or (C). The 

remaining region of the Dalitz plot (Dl) is where we expect 
to find fewer SCF events, and where the shapes for tobs and 
AE are less dependent on their position in the Dalitz plot, 
further described in Table [ml 


in the Dalitz plot of the SCF fraction and the mean differ¬ 
ence between the true and reconstructed position in the 
Dalitz plot; we include more regions in areas of the Dalitz 
plot where these quantities are largest. We use toes and 
AE PDFs specific to each region, as listed in Table lllll 
Some of the PDFs used in the parametrization of the 
SCF include Cruijff functions, Chebychev polynomials, 
Gaussian functions, and two-piece Gaussian (BGauss) 
functions. A two-piece Gaussian function is an asymmet¬ 
ric Gaussian described by the following functional form 
(omitting normalization factor) 


fBGauss{x) = exp 


— {x — to)^ 


(28) 


For the continuum background, we use an ARGUS 
function to parameterize the toes shape. The AE 
distribution is described by a linear function, and the 
BDTout distribution by an exponential function. The 
TOES: ^E, and BDTout PDFs for BB backgrounds are 
defined by the sum of the histograms from the MC simu¬ 
lations for decay modes in each background category, as 
described in Section IVl 

The continuum and BB background Dalitz plot distri¬ 
butions are included in the likelihood as two-dimensional 
histograms. For BB backgrounds, we use MC samples. 
For continuum background, we combine events from the 
off-peak data and the toes sideband in on-peak data, af¬ 
ter subtracting contributions from B decays, as described 
in Section llVl For the 2D histograms, we use the square 
Dalitz plot coordinates. A linear interpolation between 
bin centers is applied. 

The free parameters in the fit are the yields for signal. 














11 


TABLE III: List of PDFs used to describe the mss and AE 
self-crossfeed signal distributions in each of the regions of the 
^ KsTt^tv^ Dalitz plot shown in Fig. (2] The abbrevi¬ 
ations correspond to the following functional forms: Cruijff 
function described in Eq. (I27II (Cruijff), Chebychev polyno¬ 
mial (Cheb), Gaussian (Gauss), two-piece Gaussian described 
in Eq. (1281) (BGauss), and exponential (Exp). 


Dalitz plot region 

rriEs PDF 

AE PDF 


Cruijff 

Cruijff 

^kU+ (Bl) 

Cheb-b Gauss 

Exp-bSigmoid 

(B2) 

Cheb-b Gauss 

linear-bBGauss 

(B3) 

Cruijff 

Exp-bSigmoid 

(Cl) 

Cheb-b Gauss 

Cheb 

(G2) 

Cheb-b Gauss 

Cheb 

(G3) 

Cruijff 

Cheb 

(G4) 

Cruijff 

Cruijff 

Gentral region (Dl) 

Cruijff 

Cruijff 


continuum background, and BB background categories 
1 and 9. The yields for the remaining BB background 
categories are fixed to the estimated values. All the PDF 
parameters for the correctly reconstructed tues and AB 
PDFs, except for the tail parameters, are determined in 
the fit. All SCF signal PDF parameters are fixed to val¬ 
ues obtained from fits to nonresonant MC events. The 
endpoint of the ARGUS function is fixed to 5.289 GeV/c^ 
while the shape parameter is determined in the fit. The 
slope for the linear function of the AB PDF and the expo¬ 
nent for the exponential function of the BDTout PDF for 
continuum background are similarly determined in the 
fit. The isobar coefficients, x and y in Eq. (fT51) . for all 
but one of the isobar components are fitted parameters 
in the fit and are measured relative to the fixed isobar 
component. The coefficients for the reference isobar are 
fixed to a; = 1 and y = 0. In total, the fit is performed 
with 21 free parameters. 

TABLE IV: Fit fractions obtained from the fit to data when 
each additional isobar is added to the fit model one at a time. 


Additional isobar 

Fit fraction 

p(1450)+ 

0.042 

± 

0.044 

A'2*(1430)° 

0.038 

± 

0.017 

A2*(1430) + 

0.012 

± 

0.020 

A* (1680)° 

0.032 

± 

0.034 

A* (1680)+ 

0.005 

± 

0.030 


We determine a nominal signal Dalitz plot model based 
on information from previous studies (l6l Il8l - l20l| , and on 
the changes in the log likelihood in the fit to data when 
resonances are added to, or removed from, the fist shown 
in Table im In these fits to the combined B'^ and B~ data 
samples, the CP coefficients Ax and Ay are fixed to zero. 
We do not find significant contributions in the fit when 
adding the resonances p(1450)+, Ar|(1430)°, Ar|(1430)'^, 
iF*(1680)°, or Ar*(1680)+, one at a time to the default 
model. We observe that the fit fractions for these ad¬ 


ditional resonances, reported in Table HVl are consistent 
with zero. The most statistically significant of these fit 
fractions is FF (1^2 (1430)°7r+) = 0.038 ±0.017; since the 
statistical significance is less than 2 standard deviations, 
we do not include any of the additional resonances in the 
nominal fit. 

We do not observe an excess of events for invariant 
masses greater than 2 GeV/c^, suggesting that a nonreso¬ 
nant component, in addition to that included in the LASS 
parametrization, is not necessary. We observe that if we 
add a nonresonant component to the fit, the change in log 
likelihood for the binned data and the fit projections for 
the KgTT^, ATgTT^, and invariant masses are con¬ 

sistent with the expected change due to the additional 
free parameters in the fit, and do not indicate any sta¬ 
tistically significant nonresonant component. We there¬ 
fore conclude that, with the current level of statistical 
sensitivity, the base model, which includes the p{770)'^, 
Ar*(892)+, Ar*(892)°, (Ar7r)g°, and resonances, 

provides an adequate description of the data. 


VII. RESULTS 

We apply the fit described in Section |Vl] to the 31 876 
selected B'^ —>■ candidates. A first fit is per¬ 

formed on the combined B^ sample. We obtain yields of 
1014 ± 60 signal events, 24 381 ± 200 continuum events, 
2745 ± 70 BB events in category 1, and 1768 ± 140 BB 
events in category 9. The results of the fit are shown 
in Fig. [31 For the purpose of this figure, the contribu¬ 
tions of signal events are enhanced by applying the more 
restrictive selection criteria listed in Table El 

TABLE V: Selection criteria imposed to enhance the contri¬ 
butions of signal events for the results presented in Figs. O 
andjH 


Projection plot 

Selections 

mss 

-0.05 < AE < 0.05 GeV 
BDTout > 0.1 

AE 

mES > 5.27 GeV/P 
BDTout > 0.1 

BDTout 

mES >5.27 GeV/P 
-0.05 <AE< 0.05 GeV 

™-iC°7r+i ^ir+irO 

mES >5.27 GeV/c^ 

-0.05 <AE< 0.05 GeV 
BDTout > 0.1 


The branching fraction for —)> is deter¬ 

mined from the number of signal events, the efficiency 
estimated from MG events, and the total number of BB 
events in data. We take into account differences between 
the 7r° reconstruction efficiency in data and MC events, 
determined from control samples with either r leptons 
or initial-state radiation, as a function of momentum 
(^^ = 97.2%, averaged over 7r° momentum). We cor¬ 
rect for small biases in the branching fraction, as deter¬ 
mined from MC pseudo experiments generated with the 
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FIG. 3: Combined fit: Measured distributions and fit projections for candidates; (a) rriES, (b) AE, 

(c) BDTout, (d) (e) rtij^o^a, and (f) m^+^o. The points with error bars correspond to data, the solid (blue) curves 

to the total fit result, the dashed (green) curves to the total background contribution, and the dotted (red) curves to the 
continuum background component. The dash-dotted curves represent the signal contribution. The projected distributions are 
obtained from statistically precise pseudo experiments generated using the fit results. For all distributions in each panel, the 
signal-to-background ratio is increased by applying tighter selection requirement on rriES, AE, and/or BDTout, listed Table IVl 


same number of signal events and resonance composi¬ 
tion as found in the fit to data. We divide the partial 
branching fraction of 5+ —>■ if's(—>■ 7r+7r“)7r+7r° by the 
branching fraction for —>■ 7r''"7r“, and multiply the 

result by a factor of 2 to account for decay, to ob¬ 
tain the branching fraction result B (i?+ —^ itr°7r+7r°) = 
(45.9 ± 2.6 ± S.Oto'o) ^ 10“®, where the first uncertainty 
is statistical, the second is systematic, and the third is 
due to assumptions made concerning the signal model. 


The latter two uncertainties are described in Section fVlIIl 
and the breakdown of the systematic uncertainties is de¬ 
tailed in Table IxD 

We measure amplitudes and phases relative to each 
of the five two-body decays in the signal model to take 
advantage of the smaller uncertainty observed when mea¬ 
suring the relative phases of the two pairs of decays with 
same-charge K* resonances. Table I VII lists the relative 
phase, (/i, between each pair of two-body decays in the 
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TABLE VI: Combined fit: Relative phases, 4>, for the isobar amplitudes as measured from five fits to data, where each of 
the hve isobar amplitudes is in turn taken as the reference. All phases are quoted in degrees. The uncertainties are statistical 
only. 


Reference amplitude 

Resonant contribution 

A*(892)°7r+ 

Relative phase (degrees) 
A'*(892) + 7r'’ (A7r)S°7r+ (A7r);+7r‘’ 

p(770)+A° 

B+ ^ ^*(892)°^+ 


0 

-95 ± 43 

174 ± 11 

-89 ± 43 

-122 ±43 

B+ A*(892) + 7r° 


- 

0 

-90 ± 42 

6± 10 

-27 ±26 

B+ (A7r)S°7r+ 


- 

- 

0 

96 ±42 

63 ±37 

B+ 


- 

- 

- 

0 

-32 ± 25 

B+ p(770)+Ag 


- 

- 

- 

- 

0 


TABLE VII: Combined ht: Results for the fit fractions FFj (diagonal terms) and interference terms FFi^ in data for each 
resonant contribution. The uncertainties are statistical only. 



FFj and FFij 

Resonant contribution 

A:*(892)°7r+ 

A'*(892) + 7r° 

(A'7r)5°7r+ 

(A7r)o+7r° 

pinQyK'i 

B+ -> ^*(892)“^+ 

0.10 ±0.03 

0.0004 ± 0.0028 (17 ± 5) x 10”° 

0.007 ± 0.005 

-0.008 ± 0.007 

B+ A*(892)+7r° 

- 

0.14 ±0.02 

-0.010 ±0.007 

(-3± 1) X 10"® 

0.012 ±0.008 

B+ (A'7r)S°7r+ 

- 

- 

0.36 ±0.05 

(1.5 ±6.1) X 10“® 

-0.04 ±0.02 

B+ Ik-r)1+-k° 

- 

- 

- 

0.27 ±0.03 

-0.02 ±0.02 

B+ p(770) + A'° 

- 

- 

- 

- 

0.19 ±0.04 


signal model and its uncertainty. The statistical uncer¬ 
tainty in the relative phase is smallest (~ 10°) for the 
resonances that decay to the same-charge Kir state. This 
is due to a larger overlap in the Dalitz plot between the 
same-charge K* resonances than occurs for other pairs of 
resonances that only overlap in the corners of the Dalitz 
plot. 

Since the statistical uncertainties of the fit fractions do 
not depend on the reference mode, we quote in Table IVTIl 
only the fit fractions from a fit relative to the K* (892)°7r+ 
amplitude. The fit fractions for the ArQ(1430)°7r+ and 
ArQ(1430)+7r° modes are the product of the {KTr)^ S- 
wave fit fraction, shown in Table imi and the fraction 
due to the resonant contribution in the LASS parametri- 
sation (88%). The off-diagonal fit fractions are small 
compared to the diagonal elements. We calculate the 
branching fractions for the resonant contributions shown 
in Table IVIIII as the product of the total branching frac¬ 
tion and the ht fractions returned by the ht to data, 
including appropriate Clebsch-Gordan coefficients. 

To determine the overall CP asymmetry as well as the 
CP asymmetries for the contributing isobar components, 
we simultaneously ht the separate and B~ data sam¬ 
ples. The overall Acp value is calculated from the inte¬ 
grals of the positive and negative signal Dalitz plot dis¬ 
tributions. The Aa; and Ay parameters from Eq. HID are 
allowed to vary in the ht for all components except the 
reference isobar, for which the Ay parameter is hxed to 
zero (the relative phase of the 5+ and B~ Dalitz plots 
cannot be determined since they do not interfere). To 
account for possible differences in the reconstruction and 
particle identihcation efficiencies for B^ and 5“, the ef- 
hciency map as a function of the Dalitz plot position is 


TABLE VIII: Measured branching fractions B from a ht to 
the combined data sample, and CP asymmetries Acp 
(Eq. (1191) 1. The hrst uncertainty is statistical, the second is 
systematic, and the third is due to the signal model. 


Decay channel 

B(10-®) 

Acp 


45.9 ± 2.6 ± 3.0t^;g 

0.07 ± 0.05 ±0.03t2 Jl3 

A* (892)%+ 

14.6 ± 2.4 ± 1.4+°ol 

-0.12 ± 0.21 ±0.08lo °i 

A*(892)+7r° 

9.2±1.3±0.6+o;g 

-0.52 ± 0.14 ±0.04+° ®2 

Ao*(1430)%+ 

50.0 ± 4.8 ±6.lt|6 

0.14 ± 0.10 ±0.04t°+® 

Ao*(1430)+7r° 

17.2 ± 2.4 ±1.5t?;^ 

0.26 ± 0.12 ± 0.08t!j;J^ 

p(770) + A° 

9.4±1.6±1.i15;^ 

0.21 ±0.19 ±0.0712;?^ 


determined separately for B'^ and B~. The asymmetry 
for the continuum background is allowed to vary in the 
ht. The CP asymmetries of the BB backgrounds are ex¬ 
pected to be small and so are hxed to zero in the nominal 
ht. They are varied within reasonable ranges based on 
world average experimental results in order to deter¬ 
mine the associated systematic uncertainty. 

We hnd an overall CP asymmetry of 
Acp {B+ K°Tr+Tr°) = 0.07 ± 0.05 ± 0.03t^;^3, 

where the hrst uncertainty is statistical, the second is 
systematic, and the third is due to the signal model. 
This is consistent with zero CP asymmetry. Invariant 
mass projections for the ht to data allowing for direct 
CP violation are shown in Fig. 21 

Table IVIIII shows the results for the branching frac¬ 
tions and CP asymmetries obtained from the ht to data. 
The hrst uncertainty is statistical, the second is system- 
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FIG. 4: The CP fit: Measured distributions and fit projections for ^ (left column) and B K^tv tt® (right 

column) candidates; (a) mjfO^+, (b) rrij^o^-, (c) .^o(from B'^ —>■ KgTv'^Tv^), (d) mj^o^o(from B~ —^ KgTT~Tv°), (e) m^+^o, 

and (f) The points with error bars correspond to data, the solid (blue) curves to the total fit result, the dashed (green) 

curves to the total background contribution, and the dotted (red) curves to the continuum background component. The 
dash-dotted curves represent the signal contribution. The projected distributions are obtained from statistically precise pseudo 
experiments generated using the fit results. For all distributions in each panel, the signal-to-background ratio is increased by 
applying the tighter selection requirements on rriES, Al?, and/or BDTout, listed in Table lYl 


atic, and the third is the uncertainty associated with the 
signal model. We observe a significant asymmetry be¬ 
tween the and distributions in the region 

of the it'*(892)+ resonance; see Figs. lUJa) and (b). We 
determine the statistical significance, S, of a non-zero CP 
asymmetry in ^ iF*(892)+7r° from the difference be¬ 
tween the best-fit value of the likelihood, and the 


value when the CP asymmetry is fixed to zero, Cq: 

S = ^-2\n[Co/CAcp)- (29) 

Using this method, we measure a statistical signifi¬ 
cance of 3.6 standard deviations for a non-zero Acp in 
5+ —>■ iF*(892)+7r°. We obtain a consistent result of 
3.7 standard deviations for the statistical significance 
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x±A X 


FIG. 5: CP parameters {x±Ax, y±Ay) obtained from the fit 
to data for —>■ A'*^(892)7r° resonant decay including the 1 

and 2 standard deviation contours (solid and dashed curves). 
The contours are estimated by calculating the uncertainty 
and correlation between the two CP parameters. The stars 
indicate the central values of the CP parameters and the cross 
sign the origin of the plot. 


by dividing the central value of the CP asymmetry 
by the statistical uncertainty, indicating that the log- 
likelihood function is close to parabolic. Figure [S] dis¬ 
plays the contours in the complex plane of the coef¬ 
ficients c = {x + Ax,y + Ay), defined in Eq. (|1]), for 
5+ —7> iF*(892)+7r° decays, and of c = (x — Ax,y — Ay), 
defined in Eq. dS]), for B~ —7> iF*(892)“7r° decays. For 
other resonances the CP asymmetry is within 2 standard 
deviations of zero. 

We also express the complex isobar coefficients c and 
c of Eq. (fT5|) in terms of amplitudes and phases, 

c = A+C*+, (30) 

c = A_e-^^- . (31) 

Table na presents the results, measured with respect to 
the —>■ iF*(892)°7r^ reference amplitude. The statis¬ 

tical uncertainties of the separate B+ and B~ decay am¬ 
plitudes, A+ and A-, vary between 0.1 and 0.3. We thus 
obtain significant statistical precision for these terms. 
With respect to the phases, (j)+ and <))_, only the (iF7r)g° 
amplitude yields a statistically precise result. For the 
other amplitudes, the statistical uncertainty ranges be¬ 
tween 70° and 170°, and only the statistical uncertainty 
is quoted. For the more precisely determined variables, 
systematic uncertainties are evaluated as well. For the 
phases of the B^ —>• (iF7r)Q^7r° decays relative to the 
B^ —>■ iF*(892)^7r° amplitude, we obtain 

0+ {(K7v)*o+n°) - 0+ (lf*(892)+7r°) = (-14 ± 18 ± 9t^)° , 

(32) 

(P- {{Kn)l-7i°) - cf)- (A'*(892)-7r°) = (ll ± 19 ± lOtg'')” . 


VIII. SYSTEMATIC UNCERTAINTIES 

We evaluate systematic uncertainties to account for 
effects that could affect the branching fractions, phases, 
and asymmetries, by varying the fixed parameters. The 
systematic uncertainties described in this section are 
summarized in Tables IXJ through IXVIl of Appendix El 

The uncertainties associated with the branching frac¬ 
tions are listed in Table m To estimate the uncertainty 
related to the modeling of the SCF PDFs, we implement 
a simpler model consisting of only four regions in the 
Dalitz plot. The PDFs are redefined using MC events to 
match the distributions found in the newly defined re¬ 
gions. We then fit the data using the new SCF model 
and take the uncertainties to be the change in the fit pa¬ 
rameters compared to those obtained from the nominal 
fit to data. All relative systematic uncertainties due to 
the SCF Toes and AE PDFs range from approximately 
1% to 4%, except for the relative systematic uncertainty 
for the B~^ —>• p{770)'^K° decay, which is 7.5%. This is 
consistent with expectations from simulation that more 
than half the 5+ —^ p(770)+iF° events are due to SCF. 

The uncertainties associated with the number of BB 
background events are evaluated by varying the estimates 
within their uncertainties, which are primarily due to un¬ 
certainties in the branching fractions. The uncertainties 
related to the BB background toeSj AE, and BDTout 
PDFs are accounted for by varying the histogram bin 
contents according to their statistical uncertainties. The 
uncertainty is then taken as the RMS of the distribution 
of the difference in the fit parameters. The uncertain¬ 
ties related to the limited statistical precision of the MC 
and data-sideband samples are similarly accounted for by 
varying the results in the corresponding histogram bins 
by their uncertainties. 

The uncertainty in the BDTout histogram PDFs for 
correctly reconstructed and SCF signal events is deter¬ 
mined by varying the bin contents in accordance with 
the observed data/MC difference. For correctly recon¬ 
structed signal events, the tails of the asymmetric Gaus¬ 
sian PDFs for Toes and AE are fixed. To account for an 
associated uncertainty, we allow the relevant parameters 
to vary in a fit to data and use the variation in the fit 
parameters to define the uncertainty. 

To validate the fitting procedure, 500 MC pseudo ex¬ 
periments are generated, using the PDFs with parameter 
values found from the fit to data. Small fit biases are 
found for some of the fit parameters and are included in 
the systematic uncertainties. 

We also account for uncertainties in the following pa¬ 
rameters describing the signal model: the mass and width 
of each resonance and the value of the Blatt-Weisskopf 
barrier radius. The associated uncertainties are deter¬ 
mined by varying the parameters within their uncertain¬ 
ties (some of which are given in Table |lT]) and refitting. 

The uncertainties in the branching fractions related 
to particle identification, tracking efficiency, and the to¬ 
tal number of BB events are 1.0%, 1.0%, and 0.6%, re- 
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TABLE IX: Results for the relative phases (j) obtained from the combined fit, the CP amplitudes A+ and A-, and the 
CP phases </>+ and (f>- obtained from the CP fit. All parameters are measured relative to the A'*(892)°7r^ reference 

amplitude. The first uncertainty is statistical, the second is systematic, and the third is due to the signal model. Note that for 
the CP phases of all contributions except for B^ —>■ , only statistical uncertainties are quoted. 


Isobar 

4>n 

A+ 

A 

</>+ (°) 

0 

1 

A*(892)±7r" 

-95±43+^“7«o 

1.46 ± 0.22 ± 0.05+^-^| 

0.82 ±0.18 ±0.05+^-^^ 

-10 ± 112 

-98 ± 97 

(A7r)S°^± 

174 ± 11 ± 11+° 

1.74 ± 0.21 ±0.11+o?2 

2.00 ±0.27 ±0.13+0^2 

165 ± 19 ± 9^3 

190±21± iCg 

(A7r);±7r° 

-89±43li°l®7 

1.44 ± 0.22 ±0.13t2 ?J^ 

1.88 ± 0.25 ±0.14t°-22 

4± 111 

-109 ±92 

p(770)±A'“ 

-122±43t°?li« 

1.24 ± 0.01 iO.OOt^-^? 

1.54 ± 0.01 ±0.09t|^ o^ 

-50 ± 168 

-120 ±71 


spectively. We estimate systematic uncertainties in the 
branching fractions associated with the 7r° and recon¬ 
struction efficiencies to be 1.0% and 1.1%, respectively. 

Uncertainties from all the above sources are added in 
quadrature to yield the total systematic uncertainties, 
which are listed in Table IXll 

We determine changes in the branching fractions, 
AB, when the signal model is varied. The system¬ 
atic uncertainties in the branching fractions due to the 
parametrization are estimated by replacing the 
LASS model with another phenomenologically inspired 
parametrization [d^. We take the differences in branch¬ 
ing fractions with respect to the nominal fit as the sys¬ 
tematic uncertainty. This is the largest contribution to 
the uncertainty due to the model. Another uncertainty 
reflects any changes in the fit parameters for the nomi¬ 
nal model when including components that are omitted 
in the nominal fit, such as the /9(1450)''', Ar|(1430)°, and 
Ar|(1430)“''. Positive and negative variations are added 
separately in quadrature to obtain the systematic uncer¬ 
tainties due to the signal model, listed in Table IXll 
We determine systematic uncertainties in the phases 
averaged over and B~ decays from the same sources 
as considered for the branching fractions. The variations 
in the phases are measured relative to the Ar*(892)°7r+ 
amplitude. Since the differences between positive and 
negative shifts in the phases, shown in Table IXIIi are 
large in some cases, we quote for those phase shifts asym¬ 
metric systematic uncertainties. 

Reconstruction and particle identification efficiencies 
cancel to first order in the fit to CP asymmetries; there¬ 
fore the only uncertainties that are included for Aqp are 
those coming from the fit and signal model. In addition 
to this, we do not evaluate any of the uncertainties that 
are found to be negligible for the branching fractions. 

An additional uncertainty for Acp arises from having 
fixed the CP asymmetries for individual BB background 
components to the mean asymmetry averaged over all 
such components. We take the largest variation of each 
background asymmetry as the corresponding uncertainty. 

The uncertainty related to the efficiency model is deter¬ 
mined by exchanging the efficiency maps for the positive 
and negative Dalitz plots and refitting the data. We then 
take the difference in CP asymmetry with respect to the 
nominal fit as the uncertainty. 

We list in Table IXlIIl the systematic uncertainties asso¬ 
ciated with the signal CP asymmetries and the variations 


in the asymmetry due to changes in the signal composi¬ 
tion. 

We evaluate systematic uncertainties for the CP 
amplitudes and CP phases from the same sources 
as for the CP asymmetries. We list the varia¬ 
tions to the amplitudes A_|_ in Table IXIVI and to 
the amplitudes A_ in Table IXVl including the un¬ 
certainties due to changes to the signal model. Ta¬ 
ble IXVIl lists the systematic variations and model uncer¬ 
tainties for </>+(-) ^(iT7r)Q° — Ar*(892)°7r+(“)^ and 

4 >+(-) TfO - A:*(892)+(-)7r0) . 

IX. SUMMARY AND CONCLUSION 

The measured branching fractions and CP asymme¬ 
tries are summarized in Table IVIIIl and the amplitude 
and phase values in Table El including statistical, sys¬ 
tematic, and model uncertainties. We have measured for 
the first time the branching fraction and CP asymme¬ 
try for the decay B~^ —>■ A^tt+tt^. We obtain first evi¬ 
dence for direct CP violation in the intermediate decay 

Ar*(892)+7r°, with a total significance of 3.4 stan¬ 
dard deviations determined by adding statistical, system¬ 
atic, and signal-model uncertainties in quadrature and 
dividing the measured Acp by the total uncertainty. 

In addition, we have measured the branching fractions, 
CP asymmetries, and relative CP-averaged phase values 
of the decays B+ A:*(892)°7r+, B+ A:*(892)+7r°, 
B+ A:^(1430)O7r+, B+ (1430)+7rO, and 

P+ — p{770)'^K'^. The results for the branching 
fractions and CP asymmetries for P+ —>■ Ar*(892)°7r+ 
are consistent with the previous measurement from 
P+ —^ Ar+7r“7r+ decays by the Belle and BaBAR Col¬ 
laborations and the results for P+ —^ ArQ(1430)°7r+ 
are within two standard deviations from the previ¬ 
ous BaBAR measurement [E [I3- The branching 
fraction for P+ -A Ar*(892)+7r° are consistent with 
the previous measurements from the BaBAR Collab¬ 
oration for the decay mode B'^ -A and the 

result for Acp is within 2 standard deviations of the 
previous measurement [El- The branching fraction 
and Acp results for P+ —>• p(770Y+ supersede the 
previous BaBAR measurements [lj|. The CP asym¬ 
metries of B+ -A A:*(892)07r+, B+ -A (1430)°7r+, 
and B^ -A p{77Q)^K^ are all consistent with zero. 
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TABLE X: Combined measurements of branching fractions and CP asymmetries from (this analysis) and from 

BABAR analyses of B+ A'*(892)°7r+ and B+ 770(1430)%+ from B+ A'+7r"7r+ [T3|, and of B+ A'*(892)+7r° from 
B+ K+n°Tr° [3. The first uncertainty is statistical and the second is systematic. 


Decay channel 

13(10-°} Acp 

K* (892)'^ 71+ 
A'*(892) + 7r'’ 
A:o*(1430)%+ 

11.6 ± 0.5 ± 1.1 0.025 ± 0.050 ± 0.016 
8.8 ±1.0 ±0.6 -0.39 ±0.12 ±0.03 

43.4 ± 1.2 ± 5.8 0.040 ± 0.033 ± 0.033 


as expected. We obtain the first measurements 
of the branching fraction and CP asymmetry for 
B+ —7> ATq (1430)+7r°, with a significance of 5.4 standard 
deviations for the branching fraction. 

We combine our results for the branching frac¬ 
tions and CP asymmetries of B+ -A A'*(892)°7r+, 
B+ -A A:^(1430)°7r+, and B+ -A K*{892)+t:° with the 
previous BaBAR measurements. The statistical uncer¬ 
tainties and all systematic uncertainties for the CP asym¬ 
metries are uncorrelated between the measurements. For 
the branching fractions, we account for possible corre¬ 
lations when combining the systematic uncertainties. If 
the systematic uncertainties are asymmetric, the aver¬ 
age systematic uncertainty is calculated from the largest 
limit. The combined results from BaBAR for these decay 
modes are presented in Table B 

Using the world averse value for direct CP violation 
in B^ -A iT*(892)+7r“ Q and the final BaBAR result for 
direct CP violation in i3+ —>■ A'*(892)+7r°, we calculate 
AAcp for the K*7r system to be 

AAcp {K*Tr) = Acp {K*+tt°) - Aqp {K*+tv-) 

= -0.16 ±0.13. (33) 

Thus the value of AAqp in K*tt is found to be consis¬ 
tent with zero. The uncertainty in the AAcp{K*tt) result 
remains large, rendering the comparison to AAcp{Ktt), 
given in Eq. ©, inconclusive at present and motivating 
improved determinations in future experiments. 
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Appendix A: Tables of Systematic and Model 
Uncertainties 

Table m lists the uncertainties in the branching frac¬ 
tions due to systematic effects, efficiency corrections, and 
changes to the signal model. Table IXlD lists uncertainties 
in the relative phase values (for and B~ decays com¬ 
bined) due to systematic effects and changes to the signal 
model. Tables IXllIl iXIVi IXVi and IXVil list the system¬ 
atic and signal model uncertainties for the CP asymme¬ 
tries, the amplitudes for the B'^ and B~ Dalitz plots, 
A+ and A_, respectively, and the corresponding phases 
(j)+ and (j)- for the B~^ -A ATg (1430)°7r+ amplitude rela¬ 
tive to that for 5+ —>• A'*(892)°7r+, and the phase values 
for the B^ -A amplitude relative to that for 

B+ -A A:*(892)+7r°. 
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TABLE XI: Combined fit: Systematic uncertainties for the branching fraction measurements, including uncertainties due 
to the signal model. 


———___________^Resonant contribution 

Source -—___________ 

Relative Variations of branching fraction (%) 
Inclusive ^*(892)° ^*(892)+ A:o(1430)° 7^0(1430)+ 

p(770)+ 

Correctly reconstructed mss and AB PDF (Hxed parameters) 

0.8 

1.1 

0.6 

1.1 

0.7 

1.2 

Correctly reconstructed and self crossfeed signal BDTout PDFs 

3.3 

3.3 

3.4 

3.4 

4.2 

4.0 

Self crossfeed signal mss and AE PDF models 

3.0 

4.3 

3.1 

1.3 

1.8 

7.5 

Fit bias 

0.3 

0.9 

0.6 

0.5 

0.7 

0.9 

BB background mss, AE and BDTout PDFs 

0.3 

0.4 

0.2 

0.3 

0.5 

0.6 

BB background yields 

0.7 

1.2 

0.6 

0.9 

2.0 

1.8 

Background model in Dalitz plot 

1.5 

3.7 

2.8 

2.8 

2.7 

3.5 

Signal efficiency model 

0.3 

1.8 

1.0 

0.4 

0.4 

0.8 

K*{892) mass and width 

0.1 

0.7 

0.3 

0.1 

0.2 

0.1 

1 ^ 0 ( 1430 ) mass and width 

3.2 

3.8 

2.1 

8.1 

5.5 

4.0 

p(770)''" mass and width 

< 0.1 

0.2 

0.1 

0.1 

0.2 

0.3 

Blatt-Weisskopf radius 

2.3 

4.4 

2.9 

7.4 

2.9 

3.7 

Subtotal 

6.3 

9.1 

6.6 

12.0 

8.5 

11.0 

Neutral pion efficiency 


1.0 

1.0 

1.0 

1.0 

1.0 

Kg efficiency 


1.1 

1.1 

1.1 

1.1 

1.1 

Charged particle identification efficiency 


1.0 

1.0 

1.0 

1.0 

1.0 

Tracking efficiency 


1.0 

1.0 

1.0 

1.0 

1.0 

^BB 


0.6 

0.6 

0.6 

0.6 

0.6 

Total 

6.6 

9.4 

7.0 

12.2 

8.7 

11.2 

Changes due to signal model 



AB( 

lO”*^) 



{K'k)q'^/{K-kYq^ parametrization 

+8.0 

-0.3 

-0.3 

- 

- 

-1.4 

p(1450)+ 

+2.3 

+0.3 

-0.4 

+2.7 

-0.8 

-2.0 

A:2*(1430)“ and ^5(1430) + 

+1.4 

-0.3 

+0.3 

-2.6 

-0.8 

-0.3 

A:*(1680)° and ^*(1680)+ 

+1.8 

-0.1 

-0.1 

+0.6 

-1.4 

-0.2 

Total (+) 

+8.6 

+0.3 

+0.3 

+2.7 

+0.0 

+0.0 

Total (—) 

-0.0 

-0.4 

-0.5 

-2.6 

-1.8 

-2.4 
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TABLE XII: Combined fit: Systematic uncertainties due to the fit model, fixed shapes in the parametrization, and changes 
to the signal model for the relative phases (in degrees) measured relative to the A'*(892)° amplitude. 




Systematic Variations (°) 


~~~~———Resonant contribution 

Systematic ——— 

if*(892)+7r° 

ifo*(1430)°7r+ 

ifo*(1430) + 7r° 

p(770)+ifg 

Self crossfeed PDFs and mapping 

7.4 

1.4 

9.6 

6.7 

Correctly reconstructed and self crossfeed BDTout PDFs 

1.4 

0.9 

1.3 

1.3 

BB background yields 

1.7 

0.5 

1.8 

2.1 

Correctly reconstructed mss and AE PDF 

0.4 

0.2 

0.9 

1.2 

Background DP PDF 

5.7 

2.2 

4.5 

5.4 

BB background mss, AE, BDTout PDFs 

0.3 

0.2 

0.3 

0.3 

Signal efficiency model 

0.5 

0.4 

0.7 

0.7 

Fit bias 

9.3 

7.8 

7.3 

23.5 

K*{S92) mass 

1.4 

0.2 

1.2 

1.4 

if*(892) width 

0.3 

0.1 

0.1 

< 0.1 

ifo(1430) mass 

+43 

-33 

+ 5.2 

- 4.7 

+48 

-37 

+46 

-36 

ifo(1430) width 

5.2 

4.0 

5.7 

14.6 

p(770)"'' mass 

0.3 

0.1 

0.4 

1.0 

p(770)’^ width 

0.6 

0.1 

0.6 

0.2 

Blatt-Weisskopf radius 

+ 15 
-2 

+ 1.0 
- 1.4 

+ 17 
-3 

+ 0.3 

- 0.5 

Total 

+48 

-36 

+ 11 
-11 

+53 

-40 

+55 

-47 

Changes due to signal model 

/{K'k)q'^ parametrization 

-67.0 

- 

- 

-60.3 

p(1450)+ 

-18.4 

-2.8 

-11.8 

-27.4 

A:2*(1430) 

+7.8 

-3.1 

+5.5 

+11.3 

if* (1680) 

-7.8 

-4.9 

-12.9 

+11.8 

Total (+) 

+7.8 

+0.0 

+5.5 

+16.4 

Total (—) 

-69.9 

-6.5 

-17.5 

-66.2 
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TABLE XIII: Contributions to the nncertainties in the CP asymmetries for the overall and resonant isobar contributions, 
inclnding nncertainties dne to changes to the signal model. 





Systematic Variations of Acp (%) 


——Resonant contribution 

Systematic ———— 

Inclusive A'* (892)° 

7V+ A'*(892)+7r° 

A'o*(1430)°7r+ 

A'o*(1430) + 7 r ° 

p{770)+K° 

Self crossfeed PDFs and mapping 

2.0 

6.0 

1.0 

1.2 

0.9 

5.0 

Correctly reconstructed and self crossfeed BDTout PDFs 

0.7 

1.7 

1.8 

1.6 

1.1 

2.8 

BB background asymmetries 

2.5 

1.4 

1.7 

1.8 

7.5 

2.3 

Background DP PDF 

0.7 

2.7 

2.1 

2.8 

2.0 

2.5 

BB background mss, AE, BDTout PDFs 

0.2 

0.1 

0.2 

0.3 

0.5 

0.4 

Signal efficiency model 

0.2 

3.9 

2.1 

0.1 

0.9 

1.3 

Fit bias 

0.3 

1.4 

0.8 

1.9 

1.0 

1.2 

A'*(892) mass and width 

0.1 

0.2 

0.4 

0.2 

0.2 

0.2 

A'o(1430) mass and width 

1.1 

4.4 

0.5 

3.0 

2.8 

2.2 

p(770)''" mass and width 

< 0.1 

0.1 

0.1 

0.2 

0.1 

0.1 

Blatt-Weisskopf radius 

< 0.1 

0.9 

0.3 

0.8 

1.1 

0.5 

Total 

3.4 

8.1 

4.1 

4.3 

8.0 

6.9 

Changes due to signal model 

{K7ryQ^/{K7v)Q~'~ parametrization 

-0.7 

+6.2 

+2.0 

- 

- 

-8.1 

p(1450)+ 

+3.3 

+1.5 

-3.4 

-10.5 

-11.6 

-21.3 

KI{U30) 

-0.2 

+5.7 

-1.5 

-7.5 

-2.7 

+ 14.4 

A'*(1680) 

-2.2 

+6.3 

+0.5 

+4.8 

-1.9 

+ 12.3 

Total (+) 

+2.4 

+0.0 

+3.7 

+13.0 

+12.0 

+22.8 

Total (—) 

-3.3 

-10.6 

-2.0 

-4.8 

-0.0 

-19.0 


TABLE XIV: Variations in the CP amplitnde, A+, including uncertainties due to changes to the signal model. In the fits, the 
amplitudes are measured relative to the K*{892)^ amplitude. 


———__________^Resonant contribution 

Svstematic ——- 

A'*(892) + 7r° 

Variation of A+ 

(A'7r)S°7r+ (A'7r)o+7r° 

p{770)+K° 

Self crossfeed PDFs and mapping 


0.02 

0.02 

0.04 

0.02 

Correctly reconstructed and self crossfeed BDTout PDFs 

0.01 

0.03 

0.02 

0.02 

BB background asymmetries 


0.01 

0.03 

0.07 

0.02 

Background DP PDF 


0.02 

0.06 

0.04 

0.04 

BB background mss, AE, BDTout PDFs 

< 0.01 

0.01 

0.01 

< 0.01 

Signal efficiency model 


0.01 

0.01 

< 0.01 

0.01 

Fit bias 


0.01 

0.03 

0.02 

0.03 

A'*(892) mass 


0.01 

0.01 

0.01 

< 0.01 

A'*(892) width 


< 0.01 

0.01 

0.01 

< 0.1 

{Ktt)q mass 


0.02 

0.02 

0.09 

0.06 

(Kn)^ width 


0.02 

0.06 

0.01 

0.02 

p(770)+ mass 


< 0.01 

< 0.01 

< 0.01 

< 0.01 

p(770)^ width 


< 0.01 

< 0.01 

< 0.01 

< 0.01 

Blatt-Weisskopf radius 


0.02 

0.03 

0.02 

0.02 

Total 

0.05 

0.11 

0.13 

0.09 

Changes due to signal model 

{K7v)q'^/{K 7ryQ~^ parametrization 


-0.03 

— 

— 


p(1450)+ 


< 0.01 

-0.01 

-0.02 


K^{U30) 


0.05 

-0.08 

-0.08 


A'*(1680) 


-0.02 

0.07 

-0.05 


Total (+) 



+0.07 



Total (—) 


■ebh 

-0.12 
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TABLE XV: Variations in the CP amplitude, A-, including uncertainties due to changes to the signal model. In the fits, the 
amplitudes are measured relative to the K*{892)^ amplitude. 


— ——__________^Resonant contribution 

Systematic ———_ 

If*(892)-7r° 

Variation of A_ 

(If7r)5°7r“ 

p(770)-If° 

Self crossfeed PDFs and mapping 

< 0.01 

0.05 

0.04 

0.05 

Correctly reconstructed and self crossfeed BDTout PDFs 

0.02 

0.02 

0.02 

0.03 

BB background asymmetries 

0.01 

0.01 

0.06 

0.02 

Background DP PDF 

0.02 

0.05 

0.04 

0.03 

BB background tties, AE, BDTout PDFs 

< 0.01 

< 0.01 

< 0.01 

< 0.01 

Signal efficiency model 

0.02 

< 0.01 

0.02 

0.02 

Fit bias 

0.02 

< 0.01 

0.02 

0.02 

If* (892) mass 

< 0.01 

0.01 

0.01 

< 0.01 

If*(892) width 

< 0.01 

0.01 

0.01 

0.01 

{Kt:)q mass 

0.01 

0.08 

0.08 

0.04 

{Kn)Q width 

0.01 

0.07 

0.05 

0.03 

p(770)''' mass 

< 0.01 

< 0.01 

< 0.01 

< 0.01 

p(770)''' width 

< 0.01 

< 0.01 

< 0.01 

< 0.01 

Blatt-Weisskopf radius 

0.01 

0.01 

0.05 

0.01 

Total 

0.05 

0.13 

0.14 

0.09 

Changes due to signal model 

(7^7r)o'^/(i^7r)Q^ parametrization 

-0.04 

— 

— 

-0.09 

p(1450)+ 

0.04 

0.12 

0.22 

0.22 

If2*(1430) 

0.05 

0.07 

0.05 

0.03 

If*(1680) 

-0.01 

-0.02 

-0.02 

-0.04 

Total (+) 

+0.06 

+0.14 

+0.22 

+0.23 

Total (—) 

-0.04 

-0.02 

-0.06 

-0.09 


TABLE XVI: Variations in the CP phase values (j)± (in degrees) measured for the {K-k)q^-k^ amplitude relative to the 
A'*(892)°7r^ amplitude, and for the {K-k)q^tv^ amplitude relative to the A'*(892)^7r° amplitude. 


Systematic 

1 Absolute variations of CP phase values 



(If7r)SV=" - 


(If7r)«^“ 

-If*(892)±7r" 



<!>+ 

4>- 



Self crossfeed PDFs and mapping 


0.6 

4.5 

1.6 

3.1 

Correctly reconstructed and self crossfeed BDTout PDFs 

0.9 

1.6 

1.4 

2.6 

BB background asymmetries 


1.4 

1.3 

0.7 

1.0 

Background DP PDF 


2.5 

3.0 

2.0 

2.7 

BB background mss, AE, BDTout PDFs 

0.2 

0.2 

0.2 

0.4 

Signal efficiency model 


0.3 

1.4 

0.8 

1.2 

Fit bias 


2.4 

5.2 

1.7 

0.9 

If* (892) mass 


0.2 

0.2 

0.4 

0.6 

If*(892) width 


0.3 

0.3 

0.2 

0.2 

(If7r)o mass 


6.1 

5.3 

5.9 

5.5 

(If7r)o width 


4.2 

5.0 

5.2 

4.5 

p(770)’'" mass 


0.2 

0.2 

< 0.1 

0.2 

p(770)''’ width 


0.2 

0.2 

0.1 

0.2 

Blatt-Weisskopf radius 


1.0 

3.3 

1.3 

3.8 

Total 

8.5 

11.3 

8.7 

9.6 

Changes due to signal model 

p(1450)+ 


+3.7 

-0.9 

-3.2 

+9.4 

If2*(1430) 


-2.6 

+0.5 

+3.6 

-8.8 

If* (1680) 


+1.9 

+3.3 

+2.2 

+14.6 

Total (+) 


+4.2 

+0.5 

+4.2 

+17.4 

Total (—) 


-2.6 

-3.4 

-3.2 

-8.8 



